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ABSTRACT 


This paper examines the standard methods for analog differentiation. 
It concludes that the compromised derivative circuit with paralled cap- 
acitance and resistance in the feedback loop and added resistance in the 
input leg is the most satisfactory normal circuit provided that the poles 
generated by the components are suitably located. It shows furthermore 
that multipole derivative circuits can be devised from mathematical 
approximations of the type 

m (m 

pra = = Bois + Ay, tee: ay + €& 
such circuits have greater descrimination against frequencies above a 
specific level and give more accurate derivate for frequencies which 


Guewllow relative GO this CUmeert point. 
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I. INTRODUCTION 


The importance of quantitative description has been well stated by 

eretel tS liveikiae 

"I often say that when you can measure what your are 

speaking about, and express it in numbers, you know 

something about it; but when you cannot measure it, 

when you cannot express it in numbers, your knowledge 

is of a meagre and unsatisfactory kind; it may be the 

beginning.of knowledge, but you have scarcely, in 

your thoughts, advanced to the stage of science, 

whatever the matter may be.” 

iSdayeineexperimenual engineering, in System and process conturele 

we are concerned with the measurement and use of electrical signals 
derived from a wide range of varying physical states. Transducers exist 
Wiatecheeane Convery the variacvion Of measurable gquantitites invo such 
Signals. However, at the present time, with the exception of velocity 
type pick ups, there are no transducers which give outputs proportional 
bOLtuse time rable of Change of the transduced property. Likewise im 
many systems used for process control a derivative feed back loop is 
preferable to an integral feed forward. The factors of primary impor- 
tance in connection with rate mode are that by opposing all change, rate 


x 
mode has a great stabilizing effect on at ee 


The need for a clear understanding of electrical analog derivative 
circuits is apparent. On the assumption that the signal input is clean, 
standard text books on analog methodology give circuits which will pro- 
duce derivative outputs. However, the world of reality, as we know it, 
is always noisy, and consequently this theoretical ideal input is never 
achieved. Thus, the derivative circuit responding to the noisy input 
differentiates both the noise and the pertinent signal. Since in many 
natural phenomena and control processes the variable changes are 
restricted to relatively low frequency and noise is primarily generated 
or picked up at high frequency the signal to noise ratio is worse after 


differentiation. This is due to the fact that a derivative device is by 


x 
Superscript numerals refer to the Bibliography at the end of this 
report. 


on 





definition a rate of change responding device. Standard authorities 
Menem c leat iticdualie  Vermenauure wal Nee response and of the electronic 
circuitry makes this unsatisfactory situation one which the process 
control and instrumentation engineers must learn to live with, to avoid 
or to compromise. Engineering is of course the art of compromise — but 
for such to be effective its implication must be understood. The purpose 
Of this paper @egto attempt to demonstrate how analog derivative Circuats 


can be adapted to specific performance requirements. 


Oe ge 





II. STANDARD DERIVATIVE CIRCUITS 


The circuits given in Figs la-c are those generally quoted by 

| 4 
authoritative texts as the most suitable for analog differentiation. ©?” 
Hie CircuiG Shown elie IseceLrue dilterentiatora. [his eie readily 


see from the system analysis given below. 


eo(s) = -SRCe, (s) (1) 


which may be written in the time domain 


Sree, a ee (t) (2) 
This latter equation is of course a special case of the more general 


Laplace transform theorem. 


(8) «nn a os : 


The validity of Eq. (2) is not dependent upon the form of the input signal, 
fiemerciverequires that n= 1]. Thus, the circuit willy not discrimimarc 
against any frequency, but will give the derivative of the input which it 


receives. Let us consider then its behavior when the input is defined by 


e. =A sin wt (4) 
Juial 


Cont = “WRCA cos wt ©) 
We see straightaway that the sinusoidal signal is phase shifted 1/2 


and increased in amplitude by the’ factor WRC. Consider 


25 f(t) + A sin wt (6) 


and 
ae 
eee = if (t) + wA cos wt (7) 
It is clear that unless f(t) is some simple function of wt the 


relative amplitudes of the two components are not the same in the initial 





dq) 





C) 
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function and the derivative. However, in the particular case for which 


f(t) = er’ we meet the special condition that 


-e = w(t) + wA cos wt (8) 


In analog control loops the pertinent signals are generally restricted 

to the low-frequency range but this essential data is almost always mixed 
with "noise". Noise is essential high frequency. Thus the output of 
Vitomeric itrerentiavion Circuit will be the desired derivative witm che 
noise derivative superposed and by virtue of the fact that the special 
case of Eq. (8) is not usually met, the signal to noise ratio will be 
seriously reduced. Indeed, in many applications a true differentiator 


becomes a noise amplifier, e.g., when fi is zero. 


(fie eircuit of Fig. 1b is designed to solve the implicit differ— 


ential equation 
(1 -a) pe? ae (9) 


and so to obtain an approximation to the derivative. For the special 
case when (a) is unity the actual derivative results. Clearly in this 
instance there can be no advantage in this circuit since it is more 
complex than that of Fig. la and gives an identical output. It is the 
ease ara $ QO, which is of interest. The above equation is implemented on 


a computer by first integrating both sides of Eq. (9), resulting in 


i, 
iee= az +/ Zio x (10) 
0 


Or 


15 
Z = az a “igie + xX (11) 
O 


Hiewpneanscrerenunction Of this circuit is 


S 


Z 
a cle ais + 1 (12) 





or with t= ({1-a) we may write 


S 


= ee eS) = 


It is seen from Eq. (12) or (13) that when a= 1, the true deriva- 
tive results, as we previously noted. Figure 2, shows the frequency 
response in this case, (linel). For the general case a # I Mee? bigs = 
quency response is limited for circular frequencies greater than Ww 


defined by 


IL 


W == = : (14) 


This is illustrated in Fig. 2 (line 2). Thus, we readily appreciate that 
the introduction of the single pole causes input signals of frequency 


higher than Ww to pass with a constant amplification. Hence, we may 


aL 
achiiewe a COnciderable improvement iN signal nolse ratio dhea par wie ul 
circumstance by the correct positioning of the pole, i.e., by an appropriate 


choice of the parameter Ga) 


ines cransrer function for the Circuit: givenean Pig liewis 


Ey (s) y R,C,S8 ae 
=e) (RC, 4 1)(R5C,8 ae 


HOmeulemopeccial case in which Ry and Cy are Zero sand RC, #0 this 
equation is identical to Eq. (1) and the discussion of circuit l(a) is 
relevant. In the special case when either Ry Ore C,, is zero, the 
transfer function is similar to that for 1(v) and the behavior then 
follows the preceeding discussion. When these particular conditions are 
not met the transfer function contains the zero necessary to generate the 
derivative together with two arbitrarily located poles. For convenience, 


al this stage of our discussion, we shall put R.C. = R.C. = 1. dn this 


ier: 22 
case Eq. (15) may be re-written as 
Sol) RC, “a 
@;(8) (st, + 1)* 





LINE I. 





_38 LINE 2. 
o| .& 
at S 
ST +1 
_s_ 
(ST, +1) 
2 
(W;) 
0 | 
| ee LINE 3. 
LOG, W 


FIG. 2. HIGH € LOW FREQUENCY ASYMPTOTES 
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The high and low frequency asymptotes of this transfer function are given 
in Fig. 2 (line 3). The circuit has the advantage of passing, with no 


amplification, frequencies of a and attenuating all higher frequencies. 


Jig 
(for desien purposes . 1 late Gdeclred CO Opecrave Upen signal frequencies 


up to W, With a maximum error of 1 per cent, the time constant values 


O 
should be chosen such that 


al 


10W, 


C= = 
aa 1 Roos (17) 
DINcCe UneweCrron decreases approximavely as the SdWane cheule i eQguemce., 


at a frequency of vV10w wile Gieone Voll las 


difference from wW 0 


’ 
approximately O.1 wate cent. 

(iesirequency response diagram Biven in Pig. 2 showese lean scace 
fits possible by the use of multiple poles to Cause a dérivarive type 
MPU CO GtLecriminale agammey certain frequencies and yer mallvarr 
reasonable accuracy in a restricted bandwidth. Thus, the multiple pole 
system has distinct advantages in control applications. Generally 
speaking, a control signal consists of either a de or low frequency 
Signal and a superposed higher frequency noise. In most laboratories, 
the critical noise is 60cps. Thus, in a normal laboratory control 
application the circuit 1(c) finds the most common acceptance. When it 
is used it is customary to require that the 60 cycle noise passes without 
ciiaericallOonw taes location of the poles 1s then determined in accor— 


dance with the following proceedure. 


At w = 377 rad/sec, the gain is zero db. Hence 


= ae (18) 
(St + 1) Sif 


OF 


= ] (19) 
[1 = (wt )2]° + (2wt)* 





therefore 


cee Del wee oA (20) 

1 - 2(wt)2 + (wt)* + 4(wt)2 = w2 (21) 

[1 + (wt)2]* = w (22) 
ie. 

et ae Se (23) 
Hence 

a .0515 sec (24 ) 

Jw 

Thus 

w= — = 19.4 rad/sec (25) 


ad 


It can be shown that if the corner frequency is fixed at 19.4 rad/sec 

a 1 per cent error exists at 1.94 rad/sec or 0.31 eps. Thus the circuit 
discriminates against the higher frequencies at the expense of accuracy 
Piece ChmencCm Very Slow 1 requcme,) specerum: From Fig. 2 it is clear that 
with a two pole approximate derivative circuit it is preferable to make 


the poles coincident. 





IIIT. CIRCUITS BASED ON MATHEMATICAL APPROXIMATION TO THE DERIVATIVE 


>.13 Heasabilityveer Approach 


Because of the practical problems experienced in derivative circuits 
it was decided that it might be profitable to investigate the feasibility 
of employing mathematical approximations in the construction of such 


circuits. 
The derivative of f(t) with respect to +t is defined as 


(fe oa) 
labia alk 2 
fa = —— (26) 
At +0 > 


where f, is the value of f{t) at time (+t) and is is value at 
inert st). The ratio = - f,)/At approximates to the derivative 
when At is small, the smaller At the better the approximation. 
Physically, this fundamental definition can be interpreted in the 


following way: 


The time derivative of a varying voltage is equal to a constant 
noEplted by Che change in signal amplitude which occurs an a2 small 


but fixed interval of time. 


We examine then the nature of this approximation to avderivative 
When the system is applied to a signal which consists of two components — 
a smooth signal and a superposed AC signal. If the time interval between 
the sample fy and the sample f, is correctly chosen, the AC components 
will be in phase, equal in magnitude and of the same polarity. ‘Thus, 
when the signal samples are differenced the AC components will cancel 


each other out and (f, - iy) will be approximately proportional to the 


derivative of the cian component. We have implied in this approach 
that the sampling period must be a function of the frequency of the 

AC signal. The reader will notice that there is a direct comparison 
between this approach and cere! Giirerent latiem Cech les metas 
Fourier Series. It is well recognized that if we differentiate a 
Fourier Series we usually get what is not wanted in the physical problem, 


namely, a large high frequency oscillation. This, of course, is precisely 


SG: = 





what happens when we use a classic approach to obtain the derivative of 
anoisy signal. bLanczo introduced a practical approach to such a 


problem when the suggested 


[eG +5) 8-4) & i? 


as an estimate of the derivative of f(t) at the point t instead of 
the limit process which has dubious physical significance when we 


TiMienem*elaLe 2 POUTICY  perles approximaviom. 


The question we must now answer is how can this simple approach be 
Embodiewm in 4& computing circuit. Thereesappears To Bemwodictiner poscsi— 
bilities. First, we can sample the signal at discrete intervals of time 
and difference the samples. Alternatively, we can bifurcate the signal, 
delay one component relative to the other and then recombine the two 
parts through a summing network. The output of this network will be 


proportional to the derivative. 


i. seems apparent that the first approach is@iore suited Co amdieicval 
system and the second to an analog method. It is, therefore, to the 
latter that we shall direct our attention. We begin by considering how 
Uetssaporoximacion fits into an electronic network and the first pone 
Which we clearly observe is that 1/At is the overall circuit gain. 
Thus, the system does not allow complete freedom in determining At but 
rather At must be chosen to be compatible with the characteristics of 
the amplifiers which are used. We must, of course, also consider the 
Gesemaeclatacbterisulcs we desire when assecoune Tie Caircimaiee vooas 
A cursory examination would indicate that only one time delay network 
tome Gmenedmie generale the potential difierence woltenis prepen. .]emom 
to the derivative. However, if the best discrimination against noise is 


Peiired,wtnIs 15 Nol the Case. 


S.c. two Pole, Two Leg Derivative Circuit 


The analog circuit used is shown in Fig. 3. The transfer function 


Cieuiescircuigy 1s, therefore 





FIG. 3. TWO POLE, TWO LEG DERIVATIVE CIRCUIT 


Dey hee 





if 2 i 
t 
Ep (s) : ry + RL : Ry 5 + R, (28) 
e,(s) soe) ‘ 1 (2) ; 1 
r,+R, ry +R, 
By defining 

ee cm ee 

_ F = 
ak ry a RS 1 ry A RS 
_— Co% oN Re Ren Be 

an D Sa 
2 ry aE R, 2 Rs Ly + R, 


The above equation can be written as 


Sl (5) S 
>. (ae “2, (T, - tT) St. + L)(St. + 1 (29) 


ig ik Z 


if the parameters are chosen such that 


and 

am # T (5.1) 

With the exception of the gain factor a, (T, 2 t,)5 which must be 
compensated for in the summing amplifier, the analysis of this circuit 
is similar to that for two pole approximation discussed previously. The 
difference is due to the fact that the two poles must be separated and 
hengemithe desired values of the corner fmequeneies: Cangibes: she somaune 
by graphical methods using the asympotes of the frequency response gain 


CUrve . 


Accuracy dictates that the poles should be closely spaced, but, the 
(7, - a gain factor, dictates their values. This approach demonstraves 


that mathematical approximations may have advantages. 


oe ce 





53.3. The General Case 


In a paper published in the Math eepetee) W. G. Bickley has given 
a set of tables for use in numerical computation of derivatives. In 
these tables coefficients are given for insertion into the general 


Formula 


iE 
ons = pr (ADYo <= ALY, staeerewar c AY) iB & (s25 


In these formulae, m= 1, 2, 3, ete., and reference to the order of 
the derivative of y with respect to x at the point x = ra where 


vor iy on eoeee DP 


It is significant to note that in these general formulae the con- 


dition is always met. 


» a (3a) 


This general formula indicates that if we sample a signal at discrete 
times At, 2At, etc., amplify these samples by an amount proportional vo 
Bickley's coefficients and algebraically sum, the sum so derived is pro- 
pertienal! to the neh derivative. The number of time dealys (Nee) used 
in the process determines the resulting number of poles generated in the 


overall transfer function of the circuit. 


3.4. The Three Pole, Three Leg Derivative Circuit. 


From Eq. (32), the formula for the three point approximation to a 


derivative is 


ie 
y Bap MAG, + AnYo + Ags) (34) 


Figure 4 gives the analog Implementation of the above equation. 


Setting 


oe 





in C2, SS ees &o 


Cs 


FIG.4. THREE POLE, THREE LEG DERIVATIVE CIRCUIT 
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ee gee f 

i ry + Ry Ak ry + Ry 
Cerak R R 

t= 0, = = “= 
2 2 oe eae 2 

oe SET Pua Re 

= re + R,. S re, + R. 


ties Gransier function of the Gireuit becomes 


io’ S| ee eee (35) 
2a, OT, + ] ST, + 1 ST. +] 


Based on the selection of the time constant such that ae =" ite os 
suggested by Bickley's constant interval, the transfer function can be 


written as 


eas) 5 
a a tc eves | (56) 


Meovaided uhe Circuit parameter are Chosen Such thac the following teen. 


Giterons are Sacistied: 


a, - ay + a0, = 0 (37) 
6,0, - 98,0, ate 2a,0., = (38) 
5a, - 4an0i, + 5a,,0t #0 (so) 


MiesSsaveeucsau SAavist 16s therabove wares Cauactons was founds Gem 
1, -4, 3. This ratio is identical to the coefficients listed by Bickley 


POrehismuakce POlmumci tf Crenuiaud@m | Ormula. 


Thus, a three-pole approximation has been generated using the 
mathematical finite difference method. The value of three poles over 
two is obvious, they permit greater accuracy in the lower frequency 


range. The net effect of three poles can be a -12 db/octave slope 
through the higher frequency chosen to have O db gain. 


i Oi is 24 





Figure 5 shows the theoretical frequency response of the three-pole 


approximation derivative circuit given in Fig. 4. The time constants used 


we ce 


The estimated response of the double-pole circuit with equal time con- 
stants is included for comparison. The curves illustrate clearly that 
the three-pole circuit has considerably more accuracy in the lower fre- 


quency region, although both circuits pass 60 cps without gain. 


5.9. The Four-Pole, Four-Leg Derivative Circuit 


mie improvement of the low-frequency response, resulting from the 
inclusion of an additional pole leads one to consider the practicality 
Smcmeeli—-HOLe derivative newyork wshollowing the precedure used in the 


previous sections, the desired derivative is expressed in the form 
(= + + + 


Figure 6 shows the analog implementation of this equation. Setting 





| risRie, a Re 
lk ry + R, ll ry 5 RL 
a roRoC5 . RF, R,, 
= LES ST ~ —s 7 
2 tS =: R, z R, ls ee R, 
a rRC. +, R = 
S di R. > lr. Se Ro 
+ Ty + Ry + RL ry + Ry 


the transfer function of the above circuit becomes 


S55) a a,c. ; a0, ; a C., ; a) QL re 
ete) . St, + i St, + ] St_ +1 ST), + ] 


= 1 





ee LO 
(St, +1MStStIMST, +IMS t+ 1) 


S 
| (St,+1) | 
ae 
(St5+1 (StstI (ST HMST AIS Ts +1) +1) 


oot: 





| e <) lO 20 90 100 200 377 1000 
w (RAD/ SEC) 


FIG. 5 COMPARISON OF FREQUEY RESPONSE FOR THE VARIOUS 
DERIVATIVE SYSTEMS. 
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Again chosing the time constants such that T= Ani = ye nose eee 


transfer function becomes 


Enis) - 9 
Se f S27 Gas 2 Wise eS (42) 


iG! 


if the circuit parameters are chosen such that 


WB 6a, + 4a 0, - 3a, = 0 Clg) 
26a, - 19a,%, + a - lla, = 0 (444 ) 
aj, - a0, + a0, - aA, = 0 (45) 
9a,0, - 8a,a, + Ta,0, - 6a) O, #0 (46) 


The solution of the above equations yields the following relation- 


ship between the parameters 


an, = 12a, (47) 
a0, = e7a,a, (48) 


Avetis point 1 must be noted that the ratio amon the parareverce i. 
given by Eqs. (47)-(49), is not the same as the coefficients listed by 
Bickley for the four-point derivative approximation formula. Investigation 
showed that upon using the coefficients listed by Bickley, the resulting 
transter function contained an additional zero. This Zero wassor whe 

form (en + 1) thus it would cause little change in the low-frequency 


response of the circuit. 


While Eq. (42) gives a satisfactory form for a derivative approx- 


imation, the gain requirement dictated by the 6a, coefficient re- 


25 
ita 
Guires further consideration. To achieéve an exaet derivative am ties low 
PRequene ye keanee suhils "COCtmLedelu Imola. ©Cileale Uiliiby sae Liae a, 's represent 


POveMLlOMebter SeLtings, and therefore must have values of Univ ver lea 


a6. = 





oince the aa 's represent the relative gains of the four legs of the 


circuit, the gain may be written as 


a_@ 
6 ey = 1 (50) 


because 


ad, = 27 a0 (48) 


S a 


In order to discriminate against 60 cps, as with the other circuits 


discussed, a value of tT) = .005 sec. was found convenient. Hence from 
Eq. (50): 

OO 

S 


which, since 


< 
a, =e 
gives 
a, 2 900 (52) 
But a, was defined as 
Re 
2 


and therefore 


R, 2 900(r, + R,) (55) 


Hence it follows that the overall system gain must be greater than 900. 
fF eain Of this magnitude makes such 4 Circuit iampractical for operauiens 


Where economy dictates simple, unsophisticated operational amplifiers. 


The theoretical gain curve for this network is shown in Fig. 5. As 
expected, there is greater accuracy at low frequency than in the two- 


er three-polle cases previously discussed. 


2 ae 





3.6. The Four-Pole, Two-Leg Derivative Circuit. 


the impracticability of the previously described four-pole system, 
arising from the high inherent amplification requirement, caused con- 
sideration to be given to alternate approaches. The first of these is 


Shown cli Plea (Jt Ora uh tome heii We seu 


tT. = Cr ik 


i il 
a ere _ R., 
ea 1 Ba ee 
C.r.R a 

to = Te OL a= 

ee eee 2S ae 
4 ae eee 

4 er, + R, 


ie tLhe-resimtine transfer function isefi1ven. by: 


es) = [ aj, : 8,0, | oe 

Sa) ST, apr ST, qi db (St, + 1) (St), + 1) 
Tihs thenker Tuncinlonemay bes wert commas 

e (s) tS) 

= Tey 7 alt, + %,) - (4) tt) 

in (St, +1)(St, +1)(St,+1)(St, +1 

le 2 2 ) 
(55) 

if the parameters are chosen such that the following equations are 
satisfied: 

2505) - B.C, =) (56) 

TT, - Tet = O 0) 

(T, Ms t,) oy — a t ) - O (58) 


a Lae 





q, 


in 'e) 
2 2 2 2 
@ > 


Gy ~a56 


FIG. 7 FOUR POLE, TWO LEG DERIVATIVE CIRCUIT 
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As in the previous derivative approximations, we wish to choose 
the time constants of the poles such that 60 cps is passed without 
amplification. This requirement provides another equation relating the 


15 ish 1.@., 


log 377 = log 377, + log 377%, + log 377%. + log S77%), (59) 


but it clearly does not give a unique solution. The optimum approach 
appears to be to choose two of the t values and determine the other 


two from these, graphically. Figure 8 shows the proceedure adopted. 


The graphics of Fig. 8 were derived from the high and low-frequency 
asympotes of the transfer function given in Eq. (55). ‘The procedure 
was: First construct the four-equipole envelope with w=1 at zero 
db and w = 377 at zero db as the predetermined vertices (the lines 
OA and OB have slopes of +6db/octave and -18db/octave respectively). 
second, construct the double pair envelope. This was done by choosing 
point C online OA ata position corresponding to 0) = 4O. The line 
CD then drawn with the required slope of -6ab/octave, hence point D 
was determined. Next, through point C, a line with zero slope was 
drawn and through point D a line with .-12db/octave slope. The two 
undetermined poles of the four unequal pole system are then the inter- 


cepts of any line parallel to CD onthe lines CE and ED. 


This construction method is completely general as it is merely a 
graphical solution to Eqs. (57) and (59). The values obtained from the 
example given in Fig. 8 are: 


Ge =. O25 T. = 2 OlF 


WS 


0054 1), = .009 
These values satisfy all the requirements of Eqs. (57), (58), and 


(59) and yield the following transfer function for the circuit given 
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FIG. 8. METHOD OF DETERMINING POLE LOCATIONS 
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Since as 16 & POLEHLIoMeter SetCLing wand 1s thereiore diess than 
unity, for the correct magnitude of the derivative at low frequencies, 
a, must be equal to or larger than 185. This gain is much more practical 
than the 900 obtained from the previous four-pole approximation circuit. 
With the construction method outlined in Fig. 8, however, the gain 
requirement can be adjusted since it is in effect the reciprocal of 
[(t. + T)) zs (T, + t,)]. The theoretical frequency response of this net- 


work is shown in Fig. 5 for comparison with those previously discussed. 


All circuits have been examined in the laboratory, using a Donner 
Model 3500 Analog Computer, as part of this program. Our experiments 
show that the output of each circuit for various control type inputs and 
tiesirrequency response characteristics are in agreement with our thee 
rewieal predictions. The results of the frequency response vesus did 
noc vary significantly from those shown in Pig. 5. The response of the 
two-leg, four-pole derivative device, shown in Fig. /, to a range of 
inputs is given in Figs. 9-15. Results of the other circuits showed no 
marked deviation and so are not included. Figures 9-13 show the input 
Signals, the derivative of the input and the integrated derivative for 


comparison with the input signal. 


Following the method of synthesis used in the last circuit, several 
other circuits were devised using various combinations of R-C time 
delays. ‘These circuits, along with their respective Cransfer functions 
are given in Appendix A. The analysis of each circuit follows the for- 


mat previously described. 


3.{. An Integral Approximation to the Derivative. 


In Fig. 14 we show a comparison between a Fourier series, a Fejer 
sum and Lanczo's Q@ factor curves for a 12 term approximation to a 
rectangular wave. These curves indicate that the Lanczo's Q factor 
method has a clear advantage. Although we shall not present the argument 


here, we point out that the Lanczo's Q@ method is essentially the same 
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approach as the derivative approach to which we have previously referred. 
It seems reasonable, however, to consider whether or not there may be a 


possible analogy with the Fejer sum. 


If we refer to Fig. 15 then we see that the derivative is approx- 


imately given by 


Ge - f.)/2at (oi) 
Let us introduce the median PS: Then the area 

= ae 
A, Gee + is) 5 (62 ) 


and the area 


2, = (eye) = (63) 
Thus 

eA 2 (ey -t9) # 

(A,-A)  (f, - 5) a 

as | | OR = (65) 


At 


Now the area AL can be written as the difference of two integrals, i.e., 


t-At Beas 
A, -/ Goat ai f(t )dat (66 ) 
O O 
Likewise 
t +-At 
A, - | e(e)ae : f Petal (67) 
O O 
Therefore 
t eee t-2At 
Boe Be - | e(e)at Z 2 | te Che +f ml tae (68 ) 
O O O 
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Thus it follows from Eq. (68) that the derivative df/dt can be 


expressed in terms of difference in integrais. Thus 


t t-At t-2At 


a = (at )* i Ge eke = 2 | f(t jdt +/ f(t )dt (69) 
O O O 


We examine now the nature of this approximation to a derivative 
when the system is applied to a signal which consists of two components. 
A smooth signal and a superposed AC signal. We see immediately that if 
At is correctly chosen then the integrals of the AC components are 
all identical. Thus, since the algebraic sum of the gains associated 
with these integrals is zero, they will be self-compensating so far as 
AC is concerned. We observe, too, that if this mathematical expression 
is to be transformed into an electronic circuit then the input signal 
must be trifurcated and the three arms delayed relative to each other. 
At takes on two special significances. First, it is equal to the rel- 
ative time delay between adjacent arms and secondly the inverse of its 
square is the overall circuit gain. Thus, as before the system does 
not allow complete freedom in the choice of At, instead At must be 
chosen to be compatible with the characteristics of the amplifiers which 
are used. The integrating amplifiers act as very definite smoothing 
devices for cyclic noise and the delay cireuitesare very powerful peak 
or impulse filters. As in the other approximation circuits the smaller 
the At which we can use the more nearly will the output signal 


approach the true derivative level. 


3.8. Integral Difference Derivative Circuit. 


Figure 16 shows the analog network used to implement Eq. (69). The 
practicality of this circuit is limited, however, because of the tendency 
of the input integrating amplifier to saturate for any input of constant 
polarity. Analysis of the transfer function for the network given in 
Fig. 16 shows that the amplifiers one and three may be interchanged 
without effecting the function. Thus, the circuit of Fig. 1/ gives the 
same result as that of Fig. 16, but prevents saturation. The transfer 
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FIG. 16. THREE POLE DERIVATIVE CIRCUIT. 
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FIG. 17. THREE POLE CIRCUIT GIVING 1ST AND 280 
DERIVATIVES. 
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5 r, +R, 5 Yr. + R., 
it becomes 
e.(s) aa. to 2 
0 So ee f s j (71) 
S75) SRoCF (St, it 1)(St, a Gere + 1) 
if the parameters are chosen such that 
tT, = at (2) 
aja, - a0, + a0, = 19) (73) 
5a, - 4a, + 3a,d, = 0 (74) 
6a,0, - 5a,@, + 20,0, £0 Ci), 
The solution to the above equations for AQ yields the results 
an, = 20,0, (76) 
a0 = ao, (77) 


The ratio of the relative gains of the three legs is 1,-2,+l and these 
coefficients are those given by Bickley for the three point approximation 
to a second derivative. The transfer function, Eq. Cane Shows the network 


generates a second derivative, and then integrates this to obtain the first 
derivative. The analysis of the circuit is identical to that discussed pre- 
viously for three poles. The noteworth point of this approach is that both 
the first and second derivative are generated at the same time. 


= Sop = 





IV. CONCLUSIONS 


From mathematical approximations of the type 


ay (m) , 
—— =r (Ayo t+ Ay, t ce Ay) ts (32) 


mw Ssystematie family of derivative circuits can be derived. These circuits 
are at least as satisfactory as those normally adopted and in most cases 
superior for control purposes. Their superiority arises from the fact 
that they incorporate more poles than is customarily achieved. Thus their 
"cut-off frequency line" can be made appreciably steeper than the con- 
ventional -6db/octave. The slope of this line is -6(n-1)db/octave 
where n= pole number. Associated with this there is an increase of 


accuracy in derivative value for low-frequency signals. 
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